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Abstract 



It is examined whether primordial top quarks from the Big Bang could sur- 
vive under natural conditions till present days. Hybrid quark-neutron stars, 
with dense enough core to maintain the existence of top quarks, is discussed. 
In this paper only final states are investigated. Zero temperature quark-lepton 
matter is modelled as as a gas of non-interacting particles. Local charge neu- 
trality has been assumed in order to avoid volume interactions which would 
not make thermodynamical treatment possible. 
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I. INTRODUCTION 



With the claimed discovery of the top quark at FERMILAB |^ the quark family has 
become complete and masses of all quarks are known. Now we are in the position to answer 
the question of occurrance of heavy quarks, including the top, in Nature. Such a scenario 
may belong to cosmology, more exactly the physics of Big Bang or to astrophysics, namely 
superdense quark stars. They may have dense enough cores to support the existence of 
heavy quarks and our aim is to investigate the necessary conditions for the existence and 
stability of these stellar objects. 

According to the standard theories of stellar evolution, the final state of massive stars 
can be compact, dense stellar objects, neutron stars, with radius of 10 km and mass of 
0.1 — 2 Mq. Since the redshift parameter is about 0.2 on the surface of these objects one must 
use general relativity. For a non-rotating, static stellar object the Tolman-Oppenheimer- 
Volkoff [0,0 equation expresses the condition of hydrostatic equilibrium and connects nuclear 
physics, through the equation of state (EOS), with astrophysics, through gravitation. Since 
equilibrium configurations exist for any central densities, actually even for infinite central 
density as well, a quark core may support hadron matter in neutron stars. However stability 
against oscillations still must be checked. 

It is a generally held belief in the physics community that very dense nuclear matter 
(p > 5 — 6 Pnuci) actually appears in the form of quark-gluon plasma (QGP), i.e. the 
physical degrees of freedom are no longer hadron particles but quarks and gluons. 

We assume that locally the QGP matter is electrically neutral. Without this assumption 
we would have a volume self-interaction, so a non-extensive energy variable for our system 
and no thermodynamic approach could be applied. Thenceforth one could not tell either 
what different phases exist inside the star because to obtain phase boundaries one needs 
thermodynamics. 

In Section || we discuss whether primordial top quarks could survive from the Big Bang. 
Section |ITI| is about the confinement transition of the Universe. Section |IV| deals with hybrid 
stars. Section [IV C| addresses the question of stability of these stellar objects. 

We note that this paper is only a preliminary report of an on-going research work. At 
some points serious simplifying assumptions were used to make the calculations possible. 
Our main goal was to find the top peak on the mass-density plot if it exists at all. A future 
paper will rediscuss some simplifying assumptions. 



II. PRIMORDIAL TOP QUARKS FROM BIG BANG? 

Top quarks must have been present in the very early Universe in two ways. Today's ob- 
servations suggest a positive total baryon or quark number in our Universe. At temperatures 
higher than {jrit — m^c^ a substantial part of this total quark number must have been 
in t quarks. True enough, we do not know when the present baryon number was created. 
Besides an originally matter-antimatter-asymmetric Universe, which is a not too attractive 
possibility, at least three different ideas have been seriously discussed. Namely, net baryon 
number may have been generated 
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• in a non-equilibrium situation just after the GUT symmetry breaking, i.e. at cca. 
10-35 s and 10^^ GeV temperature § 

• during the electroweak transition ( t ^ 10^^° s, T ^ 1000 GeV) via anomalous pro- 
cesses through e.g. sphalerons 

• from lepton excess produced in an assumed SUSY transition, by anomalous processes 
before the electroweak transition 

Therefore the terminus ante quem for i? > seems to be 100 picosec. However, t quarks 
were present even without a quark excess, in the form of tt pairs. Such pairs are present in 
substantial quantity at T > m^c^. Both conditions hold practically at the same temperature, 
T > 175 GeV. To estimate characteristic times, let us take a radiation-dominated k = 
Friedmann universe, 

ds^ = de - R^{t){dx^ + x^rf'f]} (2.1) 
e = W= (iV7rV30)T^ (2.2) 

where N is the effective number of helicity degrees of freedom, 1 for each light boson and 
7/8 for each light fermion pair. In the GUT symmetric phase ~ 162, and just before the 
breakdown of electroweak symmetry ^ 10. The solution of Eqs.( p.l| - p.2|) is 




where the Epi Planck energy is 1.22 10^^ GeV. Hence T ^ 100 GeV at t ^ 100 picosec. 
Therefore ti pairs vanished at 100 picosec, and excess tops, if already existed then, started to 
decay into charm and up quarks. The decay process is driven by weak interaction, therefore 
at this energy difference the lifetime is definitely shorter than 100 picosec. Consequently 
top quarks practically vanished from the Universe not later than 100 picosec = 10"^'^ s after 
the Beginning. 

III. TOP QUARKS IN THE CONFINEMENT TRANSITION 

Bottom quarks start to annihilate and decay at T ^ 4.5 GeV, which corresponds to 
cca. 10 nanosec. For a review of the early history of the Universe see e.g. 0. Decay 
times are much shorter, so they vanish practically instantaneously. The same happens 
with the c quarks at 100 nanosec. The decay/annihilation begins for s quarks, too, at 
several microsecs, but that process is interrupted at 7 microsecs by a first order confinement 
transition |^ . In this transition 3-quark groups form baryons and quark pairs form mesons. 
The fine details may depend on unevaluated details of nonperturbative QCD; anyway, serious 
density differences can be expected during the transition which lasts between cca. 7 and 12 
microsecs. One possible result is the formation of dense quark blobs. The later fate of these 
dense blobs is under debates; some authors claim that the so-called symmetric quark matter 
(equal weights of u, d and s) may be stable against hadron formation 0. But the fate cannot 
be calculated without knowing the characteristic sizes of the blobs. The transition starts 
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from nucleation cores of the new phases. The evolution is uncorrelated at distances greater 
than ct , therefore at 7 microsecs the radius of the emerging blobs will be in the range of 
some kilometers. Consequently the blobs will resemble neutron stars of the order of solar 
mass. As it is well known, a neutron star is a stationary object sustained in its own gravity, 
and density inhomogenities inside can be large. Strictly speaking the central density of a 
neutron star has no upper limit; in all calculations equilibrium configurations exist for even 
infinite central density 0. (Stability against radial oscillations is another matter and will be 
studied in due course.) So, the phase transition may regenerate the heavy quarks vanishing 
earlier. (Similarly can the gravitational collapse do so, resulting in neutron or quark stars 
at much later times.) The calculation of the equilibrium states needs full general relativity 
and the equation of state of relativistic QCD plasma. 



IV. TOP QUARKS IN HYBRID STARS 

A. The Tolman-Oppenheimer-Volkoff equation of general relativity 

In general relativity theory there always exists a locally free-falling frame in which the 
laws of special relativity govern physics. For a uniform fluid at rest in this frame the T^y 
energy-momentum tensor defines what we mean by pressure and energy density. Once T^j, 
is known in a covariant form in this frame it serves as a source term in the Einstein equation 
of GTR. 

Let us restrict ourselves to static spherical objects. Then we have 4 Killing vectors in the 
space-time, of which 3 are space- like, act on 2 dimensional transitivity surface and commute 
as S'0(3); the fourth one is timelike and commutes trivially with all of the previous three 
||1CI|| . Then, up to coordinate transformations, the most general possible line element is of 
form 

ds = exp {2v{r))dt^ - exp {2X{r))dr'^ - r^dQ (4.1) 
Now, the curvature is governed via the Einstein equation 

R^k - ^gikR = -SnGTik (4.2) 

where Rik is the Ricci tensor and Tik is the energy- momentum tensor. Henceforth c = 1. For 
a fluid the latter one has the form Tik = suiUk — P{gik — UiUk). Here e is the energy density 
and P is the pressure. Both can be calculated from the equation of state. The velocity field 
of the fluid is in our case purely temporal, and of unit length. 

= (1,0,0,0) (4.3) 

For the present symmetries eq. ([4.2|) possesses three nontrivial and algebraically indepen- 
dent components, say the 00, 11 and 22 ones. However the velocity field ( |4.3| ) is a solution 
of the equation of motion, following from the automatically vanishing divergence of the left 
hand side of eq. ( |4.2| ). So now two components will suffice. We can write them as follows: 
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8nGT\ = 8nGP= -e-\- + \) + \ (4.4) 

SnGTl = 87iGP= — e~^(u// + -vl^ + + 4.5 

^ 2 2 r 2 

Now let us write, simply 

exp (-2A(r)) = 1 - 2m(r) /r (4.6) 
Then an integro-differential equation is obtained for m(r), which we write here as 

r(r-2m(r)) ' ' 

■m{r) = in Jq e{r)r'^ dr (4.8) 

and there remains a quadrature for the function z/(r), which is trivial and will be ignored 
here. In the last equation we exploited a regularity relation lim^^o rn{r) = 0. Eq. ( [4.7|) is the 
Tolman-Oppenheimer-Volkoff (TOV) equation. If the equation of state is fixed as a function 
P = P{e) or in any equivalent parametric form, then to any central density £(0) = Eq a. 
unique solution can be obtained. 



B. Phase transition 

In our environment the low density state of quarks is hadronic matter, i.e. a mixture of 
tightly bound qqq and qq clusters. It is not quite sure that this is the lowest energy state 
indeed; there are speculations that the "symmetric strange matter", i.e. an equal weight 
mixture of m, d and s in plasma state, may be energetically preferred against neutron matter, 
and then the fact that nuclear matter remains in hadrons could be explained in terms of the 
macroscopic 5* charge of such a state. Then there may be a high potential barrier between 
the local 5 = minimum of nuclear matter and the global minimum of S* = -B/3. The 
final answer is not yet known because the calculations are rather difficult for dilute quark 
plasma with QCD being there nonperturbative. Now, observe that the "potential barrier" 
argumentation may still be valid in situations where hadronic matter transforms into quarks 
in a compression process. For the primordial quarks forming a quark star between 7 and 
12 microsecs the argument is clearly irrelevant. Namely, at that time s quarks still were 
abundant; star-sized blobs of the symmetric {uds) plasma can be formed by using up the 
existing s quarks. (This is roughly the suspected scenario of "strange nugget" formation 
rT|.) Therefore it is better not to make preconceptions about the low-density crust of the 



stars under investigation. However the conditions for phase transition are clear. Hadrons 
ha are built up from quarks qi with stochiometric numbers c^: 

^aQr K (4.9) 

r 

Then phase equilibrium holds if 

Yd'^^^r = fJ-a (4.10) 

r 

Pqgp=Ph (4.11) 
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Now, we start from inside where the matter is a quark plasma with no doubts. If, going 
outward, somewhere Eqs. ( 4.1C - 4Tl| ) holds then there a phase transition happens and the 
TOV equation is to be continued with the equation of state of the new phase. 

Since now there are two different quark chemical potentials (see Subsection |1V Dj) the 
actual phase structure may be quite complicated (for a discussion see [Il2|-) However it 
seems that, if hadronic phase exists, it is dominated by neutrons; e.g. Ref. got 7/8 part 
n vs. 1/8 part p. Then for a first approximation the phase transition is 

3q ^ n 

where one quark is u, c or t, two are d, s or b; and Conds. ([4. 101 - ^^111) read as 



Hu + '^f^d ~ /in 

p{q) ^ p{n) 



(4.12) 
(4.13) 



This will be our condition to switch to nuclear equation of state. 



C. The stability of configurations 

An equilibrium configuration may still be unstable against oscillations. The sufficient 
and necessary conditions of stability against oscillations can be found in the B Appendix of 
Ref. IP . In that work a series of theorems has been proven for those who are interested only 
in the question of the stability and not in oscillation frequencies. 

The oscillation frequency squares form a series of disjoint values flj = ^"^{bq) where Eq 
is the central energy density of the configuration. A mode becomes unstable where its fl'f 
becomes zero from above. A mode can change its stability only in points where M{eo) has 
an extremum. Ignoring possible inflexion points, one gets the scheme that (in parentheses 
the shorthand terms referring to directions of turns on an R{M) diagram) for maximum 

• with R{eo) increasing (left turn on the right hand side) : the last unstable mode 
becomes stable; 

• with R{6o) decreasing (right turn on the right hand side) : the first stable mode 
becomes unstable; 

for minimum 

• with -R(£o) increasing (right turn on the left hand side) : the first stable mode becomes 
unstable; 

• with R{eQ) decreasing (left turn on the left hand side) : the last unstable mode becomes 
stable; 



(See also Ref. ||T3|-) Since until now all calculations for the neutron star peak gave stability 



on the slope upwards, stability of all peaks can be read off from the M(/i) and -R(/i) curves 
if they are sufficiently reliable in details. 
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D. The equation of state 



At first we assume that the core of the superdense quark star can be described as non- 
interacting Fermi gas of quarks and leptons at zero temperature. Having no internal energy 



production, these stellar objects will quickly cool to near absolute zero ||14 |. 

Thus the degrees of freedom taken into account are fermionic. The role of other bosonic 
degrees of freedom (e.g. quasi-particles made up of interacting fermions) will be discussed 
in a subsequent work. At zero temperature we need not deal with antiparticles either. The 
effect of nonperturbative interaction among quarks will be simulated by introducing the bag 
constant, B, which is important at lower densities to determine the interface between the 
quark core and the hadronic surface layer. This value of S = 260 MeV is derived from 
phenomenological high energy physics and may not be quite suitable for matter at these 
densities but this is the best value available now. 

Since on cosmic time scale the weak interactions will assure chemical equilibrium, one can 
introduce chemical potentials for all particle species and linear equations will hold among 
them. The chemical equilibrium is established by the following one-W exchange processes: 

s + u ^ d + u 
b + c ^ s + c 

t + b^c + b (4.14) 

These equations imply the following relations for the chemical potentials of the differently 
charged members of the quark dublets: 

IJ'u = fJ'c = fJ't = (4-15) 
fj^d = fJ^s = fJ'b = (4.16) 

However, upper heavier quarks decay into lower lighter ones via the reaction 

Q ^q + W ^q + l + ui (4.17) 

(Top quarks decay almost exclusively via t — > b+W.) The lepton conservation is insignificant 
in our calculation because the neutrinos have negligible partial pressure and they can escape 
soon after their production. Thus we keep only the chemical potential for electrons, muons 
and taus, which must be equal: fie = fi^ = fir = A*/- Moreover, the decay equation ( |4.17| ) 
couples the chemical potential of the quark doublets: 

jl = fi + l2i (4.18) 

Thus we have only two variables left. If we assume local electric charge neutrality then there 
will remain only one independent variable, say, fi, to describe the quark-lepton matter at 
zero temperature. 

Now let us summarize the necessary expressions for the density, energy density and pres- 
sure of the quark-lepton matter. For a particle of given type the thermodynamic quantities 
are the following (with h = c = 1): 
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=(#/ d'k = -ll-^^^.'l - m^fl- (4.19) 



Here the one-particle energy is Si = ^k"^ + mf , (yfj = 6 is the degeneracy for quarks {i = 
u, d, s, c, b, t) and gi = 2 is the degeneracy for leptons (i = e~ , r~). 

The total energy and pressure of the quark-lepton matter is the following: 

£ = E^e,e(/if-m2) + 54 (4.22) 
P = E^P^e{^^^-m^)-B^ (4.23) 

Q{x) is the Heaviside function, B is the phenomenological bag constant and the chemical 
potential can be determined from Eqns. ( [4.15| , ^rT6| , [4.18| ) as the function of fi and /i/. 

Plugging in Eqs. ( [4.22| , [4.23D into Eq. ( [4.7|) and using the mass continuity condition ( |4.8|) 
we still need another equation to connect our independent variables /i and fii. 

The condition of local charge neutrality yields the necessary expression to determine the 
quark chemical potential fi dependence of leptonic chemical potential nf. 

J2 n^g^e(/i2 - m^) sgn/ii = (4.24) 

i 

where = 2/3 is the electric charge for {i = u,c,t), qi = —1/3 is the charge for {i = d, s, b) 
and qi = —1 for leptons {i = e~,/i~,r~). 



E. The system of differential equations to be integrated 



We have Eqs. ( |4.7| , [4.8| , ^.24|) to be solved numerically in order to obtain an equilibrium 
configuration of a dense stellar object. 

This system of equations is governed by n of Eq.( [4.15D which determines the (central) 
energy density, i.e. the source of gravitation. The technical details of some numerical 
problems, due to the shoulder in the density profile, are discussed in Appendix ^. 

The phase transition into hadronic matter occurs where 



de 



dn 



Ph = Pqgp 

2^ld + /i« = 3/i„ + 2/ie 



H 



(4.25) 



We neglect the pressure coming from the leptons at this stage. 

From this point on we use a hadronic EOS with parameters = 195.7 and = 266.9 
until the surface of the stellar object is reached, i.e. where pn = 0. This T = EOS exhibits 
nuclear saturation and has been derived from a relativistic mean field theory incorporating 
neutral scalar and vector meson fields. The two dimensionless coupling constants of the 
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theory were matched to the binding energy and density of infinite nuclear matter in order 
to obtain results at all densities. 

Since we are discussing cold matter, the equations for the energy density and pressure 
are equivalent to knowing T^,^. This treatment is fully relativistic. 

We checked our model against a simple EOS, quadratic in the nuclear density, as well 
which yielded similar results. 

F. The shortcomings of our approach 

Let us think over what the main deficiencies in our treatment may be. 

As for the quark matter, we used a phenomenological bag model with a bag constant and 
current quark masses as parameters. The bag constant B, which parametrizes the difference 
of the exact and perturbative vacua (thus actually accounts for interactions among quarks) 
is measured only in hadron spectroscopy at low energies, so it is rather uncertain and this 
is inherited by both e an p. Similarly, the EOS was an asymptotically free quark plasma 
which is clearly a simplification. However to this point we can note that at 174 GeV the 
quark plasma is expected to be fairly perturbative. 

At the t peak we are already at energies well in the order of magnitude of the Weinberg- 
Salam mixing of electroweak interactions. While in itself this mixing may be of secondary 
importance for the EOS, a possible accompanying first order phase transition would not be 
insignificant since then a latent heat /compression energy of the order of magnitude j^^js 
would appear simulating a substantial another "bag constant". It is not clear yet if the 
Weinberg-Salam symmetry breaking transition is a first order transition or not. However no 
doubt that mw = 80 GeV and niz = 91 GeV, well below m^. 

For hadron matter we used a Walecka-type EOS whose validity range is p > 5 lO^^g cm~^. 
Also, it definetely cannot be relevant below ^ ^O^^g cm~^ where inverse /3-decay sets in. 
However, as it turns out, the shell where our hadronic EOS breaks down carries only some 
hundredths of solar mass. Anyway, our aim was not to get accurate numbers for masses 
or radii but to qualitatively investigate a central density regime not explored as far as we 
know. 

We neglected electron pressure in the hadronic matter which is significant at matter 
densities in the order of p ~ 10^^gcm~^ and could change our results for the radii. 

We shortcutted the subtle question of phase transition by saying neutrons dominate. 
However cca. 10 % protons (and electrons to ensure charge neutrality) and probably other 
hadronic particles, up to a few per cents, like A's and heavy mesons also contribute and they 
may change the phase structure. For a truly satisfactory discussion we would need a single 
EOS which contains, as limiting cases, both the hadronic and quark matter EOS. The less 
explicite (weaker first order or possibly second-order) the deconfinement phase transition is, 
the less acute this problem of super compressing and overrarifying is [^]. 

V. CONCLUSION 

Our results are summarized in a series of figures. The masses and radii of quark stars 
have been calculated at cca. 340 different central chemical potential values, and even then 
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numerical techniques were needed to smoothen the curves since at high central densities the 
inhomogeneities are substantial to cause numerical difficulties. 

Fig. |1| displays the total mass of the equilibrium configurations vs. central chemical 
potential (monotonous with central energy density). The curve starts at the point where 
the quark plasma appears in the center. 

Fig. ^ is the same but for radius. The curve is qualitatively similar to Fig. |l|. 

Fig. ^ composes the two earlier figures and depicts R{M). By means of this curve one 
can directly apply the stability criteria of Section |IV C| . For increasing central density the 
curve starts in the upper right corner, where quark matter appears, and ends up in the 
middle, where top Fermi gas could exist in the very inside of a hybrid star. 

Note that prior to this point the star is an ordinary neutron star with Walecka-type 
equation of state. (Some configurations are shown with dashed line.) There M{eo) is on the 
rise and this upward slope was found to be stable, (see e.g. Ref. Q). This turning point 
corresponds to the Landau- Oppenheimer-Volkoff point [IE 



At the first turning point (point 1) on Fig. ^i? is growing at a maximal M. Consequently 
the mode of lowest frequency becomes unstable. The next extremal point is a mirror turn 
(turn to the right on the left side), therefore the lowest frequency stable mode becomes 
unstable again. This point corresponds to the Misner-Zapolsky point []T7|. 

The form of the quark equation of state is qualitatively similar to the EOS of the Fermi 
gas of non-interacting neutrons of the early calculations and this may explain the qualita- 
tively similar structure of our R{M) curve to that of Ref. [^]. 

Beyond this point M is again increasing; we see the asymptotically damped oscillation 
of M{eo) (or M{n)) as demonstrated in Ref. 

However not too far from the M-Z point c quarks appear in the core, later they become 
relativistic and, as usual, the M^Sq) (M(yu)) curve stops to increase while reaching the 
relativistic Fermi regime. This fact puts point 3 to slightly beyond chemical potential 2 
GeV. Hence another stable mode becomes unstable, similarly as found earlier for pure non- 
interacting neutron matter. 

Fig. ^ shows that we lose still another stable mode at the next turning point (point 
4). Comparison to Fig. |1] gives that this turning point is located at /i = 5000 GeV, which 
corresponds to the appearance of h quarks beyond 4500 GeV. The minimum is extremely 
shallow and wide. 

The next turning point (point 5), located at 192 GeV, restores one mode as for stability. 
The peak is very small (0.01 Mq compared to the background) but rather expressed; it 
seems to belong to t quarks. 

Since no new degrees of freedom may appear in the equation of state beyond t quarks, we 
stopped the calculation after this turning point. Thenceforth the M^Sq) curve must display 
the well-known asymptotically damped oscillation with icreasing instability. 

Fig. 1^ shows the energy density profile at the t peak. Note the sudden drop in the 10 
cm range which coincides with the disappearance of t quarks. The jump of the first order 
phase transition is at cca. 2 km. 

Fig. ^ is the change of the two chemical potentials, as connected through charge neutral- 
ity, in a configuration at the t peak. The /i curve decreases smoothly and monotously but 
III shows a turning point at about 300 m. This surprising behaviour needs an explanation. 
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The u, d, c,s,e~ and /i~ number densities are shown in Fig. |^. (In this range < fi < 
2000 MeV, rir = 0.) Now observe that descending with fi the electron density must drop 
sunstantially at 1500 MeV since there the c quarks vanish and there remain two negative 
hght quark flavours with —1/3 e charge each to compensate for the single positively charged 
one with +2/3 e. However rig- must not drop to since Ug is slightly below riu or rid due 
to the higher s mass. Much later, at fiu + fJ'e ^ 150 MeV the s quarks vanish and then a 
great amount of electrons will be needed to balance the double charge of u over d. This way 
He- must slowly increase with decreasing fi^ after a minimum which happens to occur at 
/i = 1420 MeV. 

Ref. Ijl^ calculated quark stars up to charm, but without assuming local neutrality and 
a hadronic mantle. They did find a separate s peak, but no peak for light quarks. However 
the assumptions of the two models were sufficiently different to give such different results. 

We may conclude that, at least in the present approximation, no stable configuration 
has been found beyond the LOV point where the star is mainly a neutron star but with uds 
plasma in the centre. If our approximations are valid then no c, b or t quarks are expected 
to exist in the present Universe at least under natural conditions. 

The present paper is only a first approximation to the problem. We listed the serious 
simplifications in Section [IV b\ . Some points will be rediscussed in a future work. Note that 
e.g. at the starting point of Figs. [l| - § these calculations would need nonperturbative QCD 
results, practically unavailable now. 
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APPENDIX A: HOW TO AVOID THE NUMERICAL DIFFICULTIES AT i? = 

Since Eq. ( [4.7| ) is singular at r = one has to start integrating at a small positive f. 
At high densities it is necessary to take into account the mass inside this sphere of radius 
f. Therefore we applied a first-order expansion of Eq. ( [4.7| ). Then one can choose the 



starting point of the numerical integration at r = r where the initial (central) density (or 



the equivalent h ) drops to, say, 95 percents. This r can be expressed as r = y0.05/i^/C 
where C = 27rG{e{r = 0) + p{r = 0)){e{r = 0)/3 + p(r = 0)). 

APPENDIX B: HOW TO AVOID THE NUMERICAL DIFFICULTIES WITH THE 
SHARPLY DECREASING DENSITY PROFILE 

Since at densities where top quarks appear one has numerical difficulties with the inte- 
gration of Eqs. (|4.7| , [4.8| , [4.24| , [4.22| , [4.23| ), we resorted to various expansions. 



At infinite central energy density p = Kr is a solution to Eq. ( [4.7| ). One gets the value 

567rG 



oi K = by comparing the powers in r 
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We get a regular solution at r = if we require that p = ^f{r) is a regular function 
where /(r) is an unknown function we find by substituting into Eq. ( [4.7D . This way the 
density profile around r = is 

/^(^) = /^o - + — ^^r^^ + • • • (Bi) 

We used an expansion up to IG*'^ order that breaks down, at densities relevant to the study 
of top matter, within a few centimeters. From that radius we applied another expansion as 
follows. 

Performing perturbation calculation around the p = K/r"^ solution ( i.e. now we seek 
a solution of the form ^(1 + e(r)) where e is small ) which belongs to the infinite central 
density, one obtains an expansion at r = ro where tq is the endpoint of the validity range of 
Eq.(0): 



K 

p{r) = —(1 + air + 0A6A286alr^ + . . .) (B2) 



Here parameter ai remains free to allow to join expansion ([BT|) and expansion (|B2| ) in a 
continous way. 

We designated an arbitrary density (say, which belongs to = 25 GeV, so that the 
relation P ^ still holds ) at which we continued the numerical solution of the set of 
equations. 
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FIGURES 



Quark-neutron hybrid star at T=0 
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FIG. 1. Mass dependence of hybrid stars on the central chemical potential fi. 
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Quark-neutron hybrid star at T=0 
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FIG. 2. Radius dependence of hybrid stars on the central chemical potential ^. 
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Quark-neutron hybrid star at T=0 
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FIG. 3. Total mass vs. radius of hybrid stars. 



FIG. 4. Total mass vs. radius in the vicinity of the top peak. 
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FIG. 5. Energy density vs. radius in the model hybrid star with central density at the the top 
peak. 
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FIG. 6. Profile of the model hybrid star with central density corresponding to the top peak. 
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FIG. 7. Particle number densities vs. chemical potential /x. 



16 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9412052vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9412052vl 
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